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Abstract — A  robust  frequency-hopping  system  with  noncoherent 
detection,  iterative  turbo  decoding  and  demodulation,  and  channel 
estimation  is  presented.  The  data  modulation  is  the  spectrally 
compact  nonorthogonal  continuous-phase  frequency-shift  keying, 
which  strengthens  the  frequency-hopping  system  against  multiple- 
access  interference  and  multitone  jamming.  An  analysis  based  on 
information  theory  provides  the  optimal  values  of  the  modulation 
index  when  there  is  a  bandwidth  constraint.  The  channel  estima¬ 
tor,  which  is  derived  by  applying  the  expectation-maximization 
algorithm,  accommodates  both  frequency-selective  fading  and 
interference.  Simulation  experiments  demonstrate  the  excellent 
system  performance  against  partial-band  interference. 

I.  Introduction 

This  paper  describes  and  analyzes  a  robust  frequency- 
hopping  system  with  noncoherent  detection,  iterative  turbo  de¬ 
coding  and  demodulation,  and  channel  estimation.  The  system 
is  designed  to  be  effective  not  only  when  operating  over  the 
additive  white  Gaussian  noise  (AWGN)  and  fading  channels 
but  also  in  environments  with  multiple-access  interference  and 
multitone  jamming. 

Noncoherent  or  differentially  coherent  demodulation  has 
practical  advantages  and  is  often  necessary  because  of  the 
difficulty  of  phase  estimation  after  every  frequency  hop.  A  com¬ 
mon  choice  of  modulation  is  orthogonal  frequency-shift  keying 
(FSK).  With  FSK,  there  is  a  trade-off  between  the  bandwidth 
Bu  of  each  frequency  channel  and  the  energy  efficiency.  For 
the  AWGN  and  fading  channels,  the  energy  efficiency  can  be 
improved  by  increasing  the  alphabet  size  q  [1],  which  is  equal 
to  the  number  of  possible  transmit  frequencies  in  the  signal  set. 
The  problem  is  that  a  large  Bu  reduces  the  number  of  frequency 
channels  available  when  the  hopping  is  over  a  spectral  region 
with  fixed  bandwidth  W.  This  reduction  makes  the  system  more 
vulnerable  to  both  multiple-access  frequency-hopping  signals 
and  multitone  jamming  [2].  A  reduction  in  Bu  is  obtained  by 
using  nonorthogonal  continuous-phase  frequency-shift  keying 
(CPFSK). 

As  an  example  of  the  importance  of  Bu,  consider  multi  tone 
jamming  of  a  frequency-hopping  system  with  q- ary  CPFSK  in 
which  the  thermal  noise  is  absent  and  each  jamming  tone  has 
its  carrier  frequency  within  a  distinct  frequency  channel.  The 
uncoded  symbol-error  probability  is  approximately  [2] 

Ps=(  — )  Burb  ( p-)  ,  BuTb  <^<WTb  (1) 

\  q  J  VW  ho 


where  Zb  is  the  energy  per  bit,  Tb  is  the  bit  duration,  ItoW  is  the 
total  jamming  power,  and  Bu  is  the  uncoded  bandwidth.  This 
equation  indicates  the  significant  benefit  of  a  small  bandwidth 
in  reducing  the  effect  of  multitone  jamming. 

Robust  system  performance  is  provided  by  using  nonorthog¬ 
onal  CPFSK,  a  turbo  code,  bit  interleaved  coded  modulation 
(BICM),  iterative  decoding  and  demodulation,  and  channel 
estimation.  The  bandwidth  of  g-ary  CPFSK  decreases  with 
reductions  in  the  modulation  index  h.  Although  the  lack  of 
orthogonality  when  h  <  1  will  cause  a  performance  loss  for 
the  AWGN  and  fading  channels,  the  turbo  decoder  makes  this 
loss  minor  compared  with  the  gain  against  multiple-access  in¬ 
terference  and  multitone  jamming.  To  implement  BICM,  which 
mitigates  fading,  a  binary  encoder  is  followed  by  a  bitwise 
interleaver  prior  to  performing  the  q- ary  modulation  [3].  The 
demodulator  exchanges  information  with  both  the  turbo  decoder 
and  the  channel  estimators.  Frequency-selective  fading  changes 
the  fading  amplitude  from  hop  to  hop,  and  the  partial-band  and 
multiple-access  interference  change  the  interference  and  noise 
during  some  hop  dwell  intervals.  Consequently,  estimates  of  the 
fading  amplitude  and  the  spectral  density  of  the  interference 
and  noise  are  computed  for  a  block  size  N  that  is  smaller  than 
or  equal  to  the  number  of  symbols  in  the  hop  dwell  interval. 
The  decoder  feeds  a  priori  information  (in  the  form  of  bit 
likelihoods)  back  to  the  demodulator  and  channel  estimator, 
in  accordance  with  the  turbo  principle. 

Frequency  hopping  with  binary  orthogonal  FSK,  a  turbo 
product  code,  and  perfect  channel  information  has  been  exam¬ 
ined  in  [4].  Frequency  hopping  with  differential  ('/-ary  phase- 
shift  keying,  iterative  decoding,  and  channel  estimation  has 
been  analyzed  in  [5],  [6].  The  proposed  system  with  non¬ 
coherent,  nonorthogonal  CPFSK  has  the  following  primary 
advantages  relative  to  the  existing  systems  with  differential 
detection  or  orthogonal  modulation. 

1 .  No  extra  reference  symbol  and  no  estimation  of  the  phase 
offset  in  each  dwell  interval  are  required. 

2.  It  is  not  necessary  to  assume  that  the  phase  offset  is 
constant  throughout  a  dwell  interval. 

3.  The  channel  estimators  are  much  more  accurate  and  can 
estimate  an  arbitrary  number  of  interference  and  noise  spectral 
density  levels. 

4.  The  compact  spectrum  during  each  dwell  interval  allows 
more  frequency  channels  and,  hence,  enhances  performance 
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Fig.  1 .  Architecture  of  receiver  for  frequency-hopping  system  with  turbo  code. 
II  =  interleaver;  IT-1  =  deinterleaver. 


against  multiple-access  interference  and  multitone  jamming. 

II.  System  Model 

In  the  transmitter  of  the  proposed  system,  encoded  message 
bits  are  interleaved  and  then  placed  into  a  matrix  B  with 
m  =  log2  q  rows  and  Nd  columns.  B  is  then  transformed 
into  a  length-TV^  vector  d  with  elements  di  £  {1, 2, ...,  q}. 
Element  di  of  d  is  found  from  column  b,  of  B  through  a 
Gray  mapping.  The  vector  d  represents  the  sequence  of  tones 
that  are  frequency-translated  by  the  carrier  frequency  of  the 
frequency-hopping  waveform  during  a  signaling  interval.  The 
modulated  signal  passes  through  an  AWGN  or  fading  channel 
with  partial-band  interference. 

As  shown  in  Figure  1,  the  receiver  front-end  dehops  the 
signal.  The  dehopped  signal  passes  through  a  bank  of  q  matched 
filters,  each  of  which  is  implemented  as  a  quadrature  pair  [1], 
[2],  The  output  of  each  matched  filter  is  sampled  at  the  symbol 
rate  to  produce  a  sequence  of  complex  numbers.  Assuming  that 
symbol  synchronization  exists,  the  complex  samples  are  then 
placed  into  an  q  x  Nd  matrix  Y  whose  ith  column  represents  the 
outputs  of  the  matched  filters  corresponding  to  the  ith  received 
symbol. 

The  channel  estimator  is  derived  under  the  assumption  that 
blocks  of  N  contiguous  symbols  within  a  dwell  interval  ex¬ 
perience  the  same  fading  amplitude  and  the  same  interference 
and  noise  spectral  density,  either  or  both  of  which  could  vary 
from  one  block  to  the  next  in  an  arbitrary  manner.  If  there  are 
N  symbols  per  block,  then  there  will  be  L  =  \Nd/N]  blocks 
per  codeword.  The  matrix  Y  can  be  partitioned  according  to 
Y  =  [Yi ,  Y2 , . . . ,  YJ ,  where  the  q  x  N  submatrix  Y £  contains 
the  received  signal  vectors  corresponding  to  the  (:th  fading 
block. 


III.  Demodulator  Metrics 


The  complex  envelope  of  a  unit-energy  q- ary  CPFSK  symbol 
waveform  with  zero  initial  phase  offset  is 


Sl{t) 


1 


ej2nlht/Ta 


0  <t<Ts, 


1=1, 2,... ,q  (2) 


where  Ts  is  the  symbol  duration  and  h  is  the  modulation 
index.  Suppose  that  symbol  i  of  a  codeword  uses  unit-energy 
waveform  sdi(t),  where  the  integer  di  is  a  function  of  the 
codeword.  If  this  codeword  is  transmitted  over  a  channel  with 


fading  and  additive  Gaussian  noise,  the  received  signal  for 
symbol  i  can  be  expressed  in  complex  notation  as 


ri{t)  =  Re  aif2£~ssdi{t)eJ,27Tfct+ei 
*  =  1)2,...,  Nd 


rii(t), 


0  <t<Ts 
(3) 


where  n,(£)  is  independent,  zero-mean,  white  Gaussian  noise 
with  two-sided  power  spectral  density  NqJ2,  fc  is  the  carrier 
frequency,  di  is  the  phase,  £s  is  the  signal  energy,  and  a, 
accounts  for  the  fading  amplitude.  The  phase  0,  is  the  phase 
due  to  the  contributions  of  the  CPFSK  constraint,  the  fading, 
and  the  frequency  offset  of  the  receiver.  One  might  consider 
exploiting  the  inherent  memory  in  the  CPFSK  when  computing 
the  metric  transferred  from  the  demodulator  to  a  trellis  decoder, 
as  described  in  [7],  However,  phase  stability  over  several  sym¬ 
bols  is  necessary,  and  the  demodulator  functions  as  a  rate-one 
inner  decoder.  More  design  flexibility  exists  if  the  demodulator 
metrics  are  computed  on  a  symbol-by-symbol  basis,  and  the 
memory  in  the  turbo  code  is  exploited  rather  than  the  memory 
in  the  modulation. 

Matched-filter  k,  which  is  matched  to  sk(t),  produces  the 
output  samples 

rTs 

yk,i  =  V2  n(t)e~j2*fcts*k(t)dt  , 

Jo 

i  =  1,2, . . . ,  Nd,  k=l,2,...,q.  (4) 


The  substitution  of  (3)  into  (4)  and  the  assumption  that  each 
of  the  {s/c(i)}  has  a  spectrum  confined  to  |/|  <  fc  yields 


Vk,', 


where 


and 


ttk,. 


—  Cli  yj Ss Cp  1  H-  Tlk^i 

(5) 
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s/2  /  ni(t)e_j27r fct sk(t)dt 

(6) 
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Since  n,(f)  is  zero-mean  and  white  and  the  spectra  of  the 
{sfe(i)}  are  confined,  it  follows  that  each  nkd  is  zero-mean. 


E[nk,in*i,i\  =  NoiPi-k 

and  the  {nkd}  have  circular  symmetry: 

E[nk,initi]  =  0. 


(8) 


(9) 


Since  n,(f)  is  a  Gaussian  process,  the  real  and  imaginary 
components  of  nkyi  are  jointly  Gaussian,  and  the  set  {nk.i} 
comprises  complex-valued  jointly  Gaussian  random  variables. 

Fet  y,;  =  [yiti . . .  yqy]T  denote  the  column  vector  of  the 
matched-filter  outputs  corresponding  to  symbol  i,  and  let 
n  =  [ni,i  ■  ■  ■  nq,i]T ■  Then  given  that  the  transmitted  symbol  is 
di,  the  symbol  energy  is  £s,  the  fading  amplitude  is  a,,  the 
noise  spectral  density  is  Nm,  and  the  phase  is  di,  y,  =  y,  +n, 
where  yi  =  E[yi\di,£s,a,i,  Vo,; ,  d,} .  Equation  (5)  indicates  that 
the  kth  component  of  yi  is 

Vk,i  ~  tXis/~£se?  1  pdi  —  k. 


GO) 
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The  covariance  matrix  of  y ,  is 

R  =  E[(  yi  -  yj(; y4  -  yz)H  \  du  aif~£s,  N0i,  0t] 

=  E[  nnH]  (11) 

and  its  elements  are  given  by  (8).  It  is  convenient  to  define  the 
matrix  K  =  R/M)?  with  components 

Kki  =  Pi-k ■  (12) 


decoder  is  interleaved  and  reshaped  into  a  m  x  Nd  matrix  V 
containing  the  a  priori  information: 


Vk,i  =  log 


p{bk,i 

p(bk,i 


1|Z \zk,i) 
0|Z \zk,i) ' 


(18) 


where  conditioning  on  Z \zk  i  means  that  the  extrinsic  infor¬ 
mation  for  bit  bk.i  is  produced  without  using  zk,i-  Since  V  is 
fed  back  to  the  demodulator, 


We  can  represent  the  conditional  probability  density  function  of 
y,  given  that  the  transmitted  symbol  is  d,,  the  symbol  energy 
is  £s,  the  fading  amplitude  is  at,  the  noise  spectral  density  is 
Noi,  and  the  phase  is  Oi  as 


P(y,  \d%,  &i  \J  £$ ,  Nq  -i,  Oi ) 


ir^N^  det  K 


exp 


-^(yt-yj  K  (yi-yj 

2V0  i 


(13) 


where  K  is  independent  of  (dj,  £s,  at.  N^.  Oi). 

An  expansion  yields 

Qi  =  (y i  -  y»)ffK-1(yi  -  yj 

=  y  "K-V,  +y^K -  2Re(y,;ffK-1yJ). 

(14) 


Equations  (10)  and  (12)  indicate  that  yi  is  proportional  to  the 
dith  column  of  K  : 

y<  =  aiy/£~se^K:4i.  (15) 

Since  K  1 K  =  I,  only  the  dith  component  of  the  column 
vector  K1  yi  is  nonzero  and 

Qi  =  yiffK”Vi  +  af£s  -  2aif£~sR e(yditie~j9i).  (16) 


For  noncoherent  signals,  it  is  assumed  that  each  Oi  is  uniformly 
distributed  over  [0,  27t).  Substituting  (16)  into  (13),  expressing 
Uditi  in  polar  form,  and  integrating  over  Oi,  we  obtain  the 
probability  density  function 


p(yi\di,aif£~s,N0i)  = 


exp  l 


y<  K  Wi+q-fs 
N0i 


X  In 


det  K 

2ai\J £s  lid,  .i 


N, 


Oi 


(17) 


where  /of  )  is  the  modified  Bessel  function  of  the  first  kind 
and  order  zero.  Since  the  white  noise  n,(t)  is  independent 
from  symbol  to  symbol,  y,  with  the  density  given  by  (17)  is 
independent  of  y i,  i  ^  l. 

Let  A  and  B  denote  the  estimates  of  A  =  Nq  and  B  = 
2a\f£As,  respectively,  for  a  dwell  interval  of  N  symbols  during 
which  a,;  =  a  and  =  Nq  are  constants.  Let  bk,i  denote 
bit  k  of  symbol  i.  Let  Z  denote  the  m  x  N(i  matrix  whose 
element  zk  i  is  the  log-likelihood  ratio  for  bk,i  computed  by 
the  demodulator.  The  matrix  Z  is  reshaped  into  a  row  vector 
and  deinterleaved,  and  the  resulting  vector  z'  is  fed  into  the 
turbo  decoder.  The  extrinsic  information  v'  at  the  output  of  the 


%k,i 


log 


p(bk,i 


p(bk,i 


1  |y<,  'yfi/jvq » v<\«fc,<) 

0lyi,7|i/jvpviW,i)' 


(19) 


where  y'  =  {A,  B}.  To  calculate  (19)  for  bit  k,  first  partition 
the  set  of  symbols  V  =  { 1 , <7}  into  two  disjoint  sets  vf 
and  Tdf,  where  vf  contains  all  symbols  labelled  with  bk  =  b. 
The  extrinsic  information  can  then  be  expressed  as 

Y,de-DmP{d\yul\i/NVwi\vk,i) 

zk,i  =  log  ^ - - - JJ, - 7 - r - y  (20) 

T,d£V(0)  P(d \yl7  7ri/jv] .  Vi\Vk,i) 

From  Bayes’  rule. 


p(d\y,j',v\vk) 


P(y\d,  y,v\vk)p(d,  Y ,  v\ufc) 
p(y,7',v\vk) 


(21) 


After  conditioning  on  d  and  7',  y  is  independent  of  v  and  thus 
p(yjd,  Y,  v\vk)  =  p(y\d,Y).  From  the  definition  of  condi¬ 
tional  probability,  p(d,  y,v\vk)  =  p(d\Y,  v\vk)p(Y,  v\ufc). 
After  conditioning  on  v\vk,  d  is  independent  of  7',  and 
thus  p(d\y,  v\vk)  =  p(rf|v\ufe).  Gathering  all  these  factors, 
inserting  them  into  (20),  and  cancelling  common  factors  yields 


Zk,%  =  log 


E  d&yy  p(yi\d’  7 \i/N^  )p(d\^i W,») 
E  deX)(<»  p(.yi\d,Yli/N]  )p{d\vi\vk,i) ' 


(22) 


This  expression  clearly  delineates  the  contribution  of  the  chan¬ 
nel  observation  y  and  channel  estimate  7',  which  influence 
only  the  p(y\d,  7')  term,  and  the  contribution  of  the  a  priori 
information  passed  to  the  demapper  from  the  decoder,  which 
affects  only  the  p(d\v\vk)  term. 

The  term  p(d\v\vk)  in  (22)  is  only  computed  for  those  d  £ 
in  which  case  p{bk  =  6|v\ufc)  =  p(bk  =  b)  =  1.  Thus, 
under  the  assumption  of  independent  code  bits  (achieved  by 
proper  interleaving),  the  probability  of  d  given  the  a  priori 
input  v  is 


p(d\v\vk)  =  X\p{bj{d)\vj),  d  G  V(kb)  (23) 
j=i 

where  bj(d)  is  the  value  of  the  jth  bit  in  the  labelling  of  symbol 
d.  The  a  priori  input  is  interpreted  by  the  demapper  to  be  v  = 
l°g[p/(l  ~  p)],  where  p  is  the  decoder’s  most  recent  estimate 
of  the  probability  that  the  corresponding  code  bit  is  a  one. 
Inverting  the  logarithm  and  solving  for  p  yields  p  =  ev /( 1  + 
ev ),  which  the  demapper  uses  for  p(b  =  l|u).  Similarly,  the 
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demapper  uses  1  —  p  =  1/(1  +  ev)  for  p(b  =  0|v).  Since 
b  =  {0, 1},  we  obtain 

m  ~Vjbj(d) 

p(d\v\vk)  =  n  !  _|_  eVj  ,  dGVk]-  (24) 

i=i 

Substituting  (17)  and  (24)  into  (22)  and  cancelling  common 
factors,  we  obtain 

E*^1)  /o  h\i/N-\  \ydui\)  nr=i  exP  (bj(d)vj:i) 

_  1  k 

~k'1  °6  J2devm  i0  (i\i/N]  \ydiA)  nr=i  exp 
k 

(25) 

where  only  the  ratio  7  =  B/A  is  needed  rather  than  the 
individual  estimates. 

IV.  Channel  Estimators 

Since  under  block  fading  and  time-varying  interference,  A 
and  B  can  change  on  a  block-by-block  basis,  each  block  is 
processed  separately  and  in  an  identical  fashion.  To  maintain 
robustness,  the  estimators  make  no  assumptions  regarding  the 
distribution  of  the  quantities  to  be  estimated,  nor  do  they  make 
any  assumptions  regarding  the  correlation  from  block  to  block. 
The  estimators  directly  use  the  channel  observation  for  the  £th 
block,  Y(,  while  the  observations  of  the  other  blocks  are  used 
indirectly  through  feedback  of  extrinsic  information  from  the 
decoder.  Since  the  form  of  the  estimation  algorithm  is  the  same 
for  each  block,  in  the  following  discussion  we  assume  that 
£  =  1  and  without  ambiguity  omit  the  dependence  on  t.  Thus  in 
this  section,  Y  is  a  generic  qxN  received  block,  d  =  \d\,  ...djv] 
is  the  corresponding  set  of  transmitted  symbols,  and  {A,  B}  is 
the  corresponding  set  of  channel  estimators. 

Rather  than  solving  the  ML  problem  directly,  the 
expectation-maximization  (EM)  algorithm  can  be  used  as  an 
iterative  approach  to  solving  this  problem  [9].  The  EM  algo¬ 
rithm  consists  of  two  steps  per  iteration,  the  expectation  step 
(E-step)  and  the  maximization  step  (M-step).  Let  {  Y,  d}  denote 
the  complete  data  set.  The  log-likelihood  of  the  complete  data 
set  is 


L(A,  B)  =  logp(Y ,  d| A,  B) 

=  log  p(Y|d,  A,  B)  +  logp(d) 

~log  p(Y\d,A,B).  (26) 


In  the  last  step,  the  logp(d)  have  been  removed  from  the 
likelihood  because  it  is  independent  of  A  and  B  and  therefore 
does  not  affect  the  maximization. 

Since  y,  and  y 1  are  independent  for  i  A  k  (17)  implies  that 


p(Y|d,  A,  B) 


exp 


D 

A 


where 


^  +  Ef=ilog/o(%^) 

(niAi  det  K)w 

N 

d  =  J2  y  "k-1*. 

i=l 


(27) 


(28) 


After  dropping  irrelevant  constants,  (27)  and  (28)  yield 


L(A,B) 


—qN  log  A  —  — 


NB2 
4  A 


N 

+ ^2  log  /() 

i= 1 


(29) 


The  form  of  this  equation  indicates  that  the  parameters  A  and 
B  must  both  be  estimated  rather  than  just  the  ratio  B/A. 

Let  r  denote  the  EM  iteration  number,  and  A^ ,  B^  the 
estimates  of  A ,  B  during  the  rth  iteration.  Applying  the  E-  and 
M-steps  with  (29)  in  a  similar  manner  to  that  in  [10],  it  is  found 


that 


A(r> 


1 

qN 


N(B M)2 
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is  the  normalization  factor  forcing  Y2k=i  Pk  i  = 
- _ 1 _ 


Y.l,Iio(^SA)p(dl  =  k) 


(33) 


and  p(di  =  k)  is  the  probability  that  dj  =  k  estimated  by  the 
decoder. 

While  a  closed  form  solution  to  (31)  is  difficult  to  obtain,  it 
can  be  found  recursively  [11],  The  recursion  involves  initially 
replacing  B ^  on  the  right-hand  side  of  (31)  with  B<'r~1') 
from  the  previous  EM  iteration.  To  select  an  initial  estimate 
for  B,  consider  what  happens  in  the  absence  of  noise.  Without 
noise,  (5)  implies  that  either  \yk,i\  =  a\jEs  (when  k  =  di)  or 
\yk,i\  =  0  (otherwise).  Thus,  an  estimate  for  ayfE~s  =  B/2  can 
be  achieved  by  taking  the  maximum  \yk,i\  over  any  column  of 
Y.  To  account  for  the  possibility  of  noise,  the  average  can  be 
taken  across  all  columns  in  the  block,  resulting  in 


B(  °) 


2  .  ' 


(34) 


The  initial  estimate  of  A  is  found  from  B>>)]  by  evaluating  (30) 
for  r  =  0.  After  the  initial  values  /1'0;  and  lk>,!  are  calculated, 
the  initial  probabilities  \ }  are  calculated  from  (32)  and  (33). 
The  EM  algorithm  terminates  when  B(r>  converges  to  some 
fixed  value,  typically  in  fewer  than  10  EM  iterations. 


V.  Selection  of  Modulation  Index 

Let  Bmax  denote  the  maximum  bandwidth  of  the  CPFSK 
modulation  such  that  the  hopping  band  accommodates  enough 
frequency  channels  to  ensure  adequate  performance  against 
multiple-access  interference  and  multitone  jamming.  We  seek 
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to  determine  the  values  of  h,  q.  and  code-rate  R  of  the  turbo 
code  that  provide  a  good  performance  over  the  fading  and 
AWGN  channels  in  the  presence  of  partial-band  interference. 
For  specific  values  of  the  modulation  parameters  h  and  q, 
the  code  rate  is  limited  by  the  bandwidth  requirement.  Let 
BuTb  denote  the  normalized,  uncoded,  99-percent  power  band¬ 
width  of  the  CPFSK  modulation.  This  value  can  be  found  for 
nonorthogonal  CPFSK  by  numerically  integrating  the  power- 
spectrum  equations  [1]  and  will  be  valid  for  frequency-hopping 
signals  provided  that  the  number  of  symbols  per  dwell  interval 
is  large  [2],  When  a  code  of  rate  R  is  used,  the  bandwidth 
becomes  Bc  =  Bu/R.  Since  Bc  <  Bmax  is  required,  the 
minimum  code  rate  that  achieves  the  bandwidth  constraint  is 
Bmin  —  BuJBrnax. 

Guidance  in  the  selection  of  the  best  values  of  h,  q,  and  R  > 
Rmin  is  provided  by  information  theory  [12].  For  specific  values 
of  h  and  q,  we  evaluate  the  capacity  C( 7)  as  a  function  of 
7  =  £S/N0  under  a  bandwidth  constraint  for  both  the  Rayleigh 
and  AWGN  channels.  Since  the  noncoherent  demodulator  will 
include  channel  estimation,  perfect  channel-state  information 
is  assumed.  Symbols  are  drawn  from  the  signal  set  with  equal 
probability.  With  these  assumptions,  a  change  of  variables  with 
u  =  yt/ \/£l,  and  (17),  the  capacity  for  the  fading  channel  may 
be  expressed  as  [3] 


C( 7)  =  log2  q  -  -  /  /  POMuk  a) 

®  v=l 


J2  J0  (207  M) 


x  log; 
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where  p(a)  is  the  density  of  the  fading  amplitude,  the  (2 q  + 
1)— fold  integration  is  over  all  values  of  a  and  the  2 q  real  and 
imaginary  components  of  u,  and 


p(u|i/,a)  = 


79  exp  [—7(11^  K  1u  + 
7 det  K 


-J0  (2a7  |u„|) .  (36) 


Equation  (35)  is  numerically  integrated  by  the  Monte  Carlo 
method.  To  determine  the  minimum  £b/Na  necessary  to  main¬ 
tain  C( 7)  above  the  code  rate  R,  we  use  the  relationship 
£s  =  R£b  log2  q  and  solve  the  equation 


R  =  C  (R£b  log2  qj  Na)  (37) 


for  all  code  rates  such  that  Rmin  <  R  <  1.  For  noncoherent 
systems  under  severe  bandwidth  constraints,  the  R  that  min¬ 
imizes  £b/Na  will  typically  bc  R  =  Rm in,  but  under  loose 
bandwidth  constraints  the  R  that  minimizes  £b/N0  could  pos¬ 
sibly  be  larger  than  Rmin  (in  which  case  the  actual  bandwidth 
is  less  than  Bmax). 

Figures  2  and  3  show  plots  of  the  minimum  £b/Na  versus 
h  for  2  <  q  <  32,  =  2,  and  BmaxTb  =  00.  When 

BmSiXTb  =  2,  the  curves  are  truncated  because  there  is  a 
maximum  value  of  h  beyond  which  no  code  exists  that  satisfies 
the  bandwidth  constraint.  Figure  2  is  for  the  AWGN  channel, 
and  Figure  3  is  for  the  Rayleigh  fading  channel.  For  each  of  the 
five  values  of  q  and  two  channel  types,  capacity  curves  were 


Fig.  2.  Minimum  £b/N0  versus  h  for  the  AWGN  channel,  2  <  q  <  32, 
BmaxL  =  2,  and  B  max  Tb  =  00. 


Fig.  3.  Minimum  £b/N0  versus  h  for  the  Rayleigh  channel,  2  <  q  <  32, 
BmaxL  =  2,  and  B  max  Tb  =  00. 


generated  for  values  of  h  ranging  from  h  =  0.01  to  h  =  1  in 
increments  of  0.01.  Thus  a  total  of  1,  000  capacity  curves  were 
generated  and  each  curve  was  created  using  at  least  2  million 
simulated  symbols  per  SNR  point  in  the  range  of  interest. 
Altogether,  nearly  1  trillion  symbols  were  simulated,  and  it 
is  estimated  that  this  task  would  have  taken  about  one  year 
to  complete  on  a  single  PC  computer.  To  speed  the  run  time, 
simulations  were  executed  on  a  virtual  private  grid  computer 
powered  by  the  idle  capacity  of  30  workstations  located  in  the 
teaching  laboratories  at  the  last  author’s  institution1.  The  entire 
simulation  scenario  took  just  two  weeks  to  complete  on  the 
grid  computer. 

'Job  scheduling  was  performed  online  via  the  Global  Grid  Exchange 
(g2ex.com)  which  runs  the  Frontier  Grid  Platform  developed  by  Parabon 
Computation  (parabon.com). 
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For  each  value  of  q,  in  each  figure  there  is  an  optimal  value 
of  h  that  gives  the  smallest  value  of  the  minimum  £b/Nn. 
This  smallest  value  decreases  with  q ,  but  there  are  diminishing 
returns  and  the  implementation  complexity  increases  rapidly 
for  q  >  8.  Let  fe  denote  the  offset  in  the  estimated  carrier 
frequency  at  the  receiver  due  to  the  Doppler  shift  and  the 
frequency-synthesizer  inaccuracy.  The  separation  between  adja¬ 
cent  frequencies  in  a  CPFSK  symbol  is  hfb / R  log2  q,  where  fb 
denotes  the  information-bit  rate.  Since  this  separation  must  be 
much  larger  than  fe  if  the  latter  is  to  be  negligible  as  assumed 
in  (4), 


fe  « 


hfb 

R  log2 1 


(38) 


is  required.  Since  the  optimal  h  decreases  while  R  log2  q 
increases  with  q,  (38)  is  another  reason  to  prefer  choose  q  <  8. 
For  q  =  4  in  Figure  3,  h  =  0.16  is  the  approximate  optimal 
value  when  BmaxTb  =  2,  and  the  corresponding  code  rate 
is  approximately  R  =  16/27.  For  q  =  8,  h  =  0.32  is 
the  approximate  optimal  value  when  BmaxTb  =  2,  and  the 
corresponding  code  rate  is  approximately  R  =  8/15.  For  both 
q  =  8  and  q  =  4,  (38)  is  satisfied  if  fe  «  0.2//.  At  the 
optimal  values  of  h,  the  plots  indicate  that  the  loss  is  less  than 
1  dB  for  the  AWGN  channel  and  less  than  2  dB  for  the  Rayleigh 
channel  relative  to  what  could  be  attained  with  the  same  value 
of  q,  h  =  1  (orthogonal  CPFSK),  and  an  unlimited  bandwidth. 


VI.  Simulation  Results 

Simulation  experiments  were  conducted  to  assess  the  perfor¬ 
mance  of  frequency-hopping  systems  with  quaternary  CPFSK 
(q  =  4)  and  octal  CPFSK  {q  =  8)  under  the  bandwidth 
constraint  BmaxTb  =  2.  The  approximate  optimal  values  of 
h  and  R  determined  from  the  bandwidth  constraint  and  infor¬ 
mation  theory  are  used.  The  interference  signal  is  modeled  as 
partial-band  noise  interference  that  introduces  Ito/ R  additional 
interference  spectral  density  in  an  interfered  frequency  channel, 
where  /i  is  the  fraction  of  the  hopping  band  with  interference 
and  Ito  is  the  spectral  density  when  //  =  1 .  Thus,  the  total 
interference  power  is  conserved  as  //  varies.  The  simulated 
system  uses  the  turbo  code  from  the  UMTS  specification  [13] 
with  2048  information  bits  and  the  specified  code-rate  matching 
algorithm.  The  receiver  executes  no  more  than  20  iterations,  as 
an  early  halting  routine  stops  the  iterations  once  the  data  is 
correctly  decoded.  The  figures  display  the  minimum  value  of 
£b/N0  necessary  to  obtain  a  bit  error  probability  equal  to  1 0  3 
versus  ft  for  Rayleigh  fading,  Ricean  fading  with  factor  K  =  10 
dB,  and  the  AWGN  channel.  A  block  coincides  with  a  dwell 
interval,  and  the  parameter  A  represents  the  spectral  density  due 
to  the  noise  and  the  interference  during  a  dwell  interval.  The 
symbols  of  a  dwell  interval  undergo  the  same  fading  amplitude, 
and  the  fading  amplitudes  are  independent  from  block  to  block, 
which  models  the  frequency- selective  fading  that  varies  after 
each  hop.  The  bandwidth  is  assumed  to  be  sufficiently  small 
that  the  fading  is  flat  within  each  frequency  channel. 

Figures  4  and  5  plot  the  results  for  quaternary  CPFSK 
and  octal  CPFSK,  respectively,  when  there  are  32  hops  per 


Fig.  4.  Minimum  required  £b/N0  of  frequency-hopping  system  with  qua¬ 
ternary  CPFSK.  h  =  0.46,  32  hops  per  codeword,  partial-band  interference, 
£b/Ito  =  10  dB,  and  bit  error  probability  =  10— 3. 


Fig.  5.  Minimum  required  £b/N0  of  frequency-hopping  system  with  octal 
CPFSK,  h  =  0.32,  32  hops  per  codeword,  partial-band  interference,  £b/It 0  = 
10  dB,  and  bit  error  probability  =  10— 3. 


codeword  and  64  information  bits  per  hop,  and  £b/Ito  =  10  dB. 
For  quaternary  CPFSK,  there  are  3456  code  bits  in  a  codeword, 
108  code  bits  per  hop,  and  54  code  symbols  per  hop.  For  octal 
CPFSK,  there  are  3840  code  bits  in  a  codeword,  120  code  bits 
per  hop,  and  40  code  symbols  per  hop.  For  quaternary  CPFSK, 
h  =  0.46  and  R  =  16/27,  whereas  for  octal  CPFSK,  h  =  0.32 
and  R  =  8/15.  Comparison  of  these  two  figures  indicates  that 
octal  CPFSK  has  a  nearly  2  dB  advantage  in  £t, /N0  relative  to 
quaternary  CPFSK  for  Ricean  fading  and  AWGN,  and  much 
more  for  Rayleigh  fading.  Both  figures  indicate  that  /r  =  1  or 
interference  over  the  entire  hopping  band  is  the  worst  case  for 
Rayleigh  and  Ricean  fading.  For  the  AWGN  channel,  a  smaller 
value  of  //  is  worst.  The  use  of  the  EM  channel  estimators  is 
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Fig.  6.  Minimum  required  £b/N0  of  frequency-hopping  system  with 
quaternary  CPFSK,  h  =  0.46,  Rayleigh  fading,  partial-band  interference, 
Sb/It 0  —  13  dB,  and  bit  error  probability  =  10-3.  N  is  the  number  of 
code  symbols  per  hop. 

shown  to  produce  a  negligible  loss  relative  to  perfect  channel- 
state  information  (CSI). 

If  the  hop  rate  increases,  the  increase  in  the  number  of 
independently  fading  dwell  intervals  per  codeword  implies  that 
more  diversity  is  available  in  the  processing  of  a  codeword. 
However,  the  shortening  of  the  dwell  interval  makes  the  channel 
estimation  less  reliable  by  providing  the  estimator  with  fewer 
samples.  Figures  6  and  7  show  the  results  for  quaternary 
CPFSK  and  octal  CPFSK,  respectively,  when  the  hop  rate  is 
varied  so  that  there  are  16,  32,  or  64  hops  per  codeword. 
Independent  Rayleigh  fading  occurs  during  each  dwell  interval, 
£b/ho  =  13  dB,  and  the  information-bit  rate  is  maintained. 
For  quaternary  CPFSK,  there  are  N  =  108,  54,  or  27  code 
symbols  per  hop.  For  octal  CPFSK,  there  are  N  -  80,  40,  or 
20  code  symbols  per  hop.  For  quaternary  CPFSK,  h  =  0.46 
and  R  =  16/27,  whereas  for  octal  CPFSK,  h  =  0.32  and 
R  =  8/15.  Comparison  of  these  two  figures  indicates  that  for 
Ricean  fading  and  AWGN,  octal  CPFSK  maintains  its  nearly  2 
dB  advantage  in  £b/Na  relative  to  quaternary  CPFSK.  Despite 
the  slow  decline  in  the  accuracy  of  the  EM  channel  estimates, 
the  diversity  improvement  is  sufficient  to  produce  an  improved 
performance  as  the  number  of  code  symbols  per  hop  decreases. 
However,  decreasing  to  fewer  than  20  code  symbols  per  hop 
will  begin  to  broaden  the  spectrum  significantly  [2]  unless  the 
parameter  values  are  changed. 

Simulation  experiments  have  confirmed  the  major  perfor¬ 
mance  advantage  of  the  proposed  frequency-hopping  system 
against  multiple-access  interference.  However,  because  of  space 
limitations,  an  explicit  demonstration  is  deferred  to  a  subse¬ 
quent  paper. 

VII.  Conclusions 

A  noncoherent  frequency-hopping  system  with  nonorthog- 
onal  CPFSK  has  been  designed  to  be  highly  robust  in  en¬ 


Fig.  7.  Minimum  required  £i/N0  of  frequency-hopping  system  with  octal 
CPFSK,  h  =  0.32,  Rayleigh  fading,  partial-band  interference,  Sb/Ito  =  13 
dB.  and  bit  error  probability  =  10— 3.  N  is  the  number  of  code  symbols  per 
hop. 

vironments  including  frequency-selective  fading,  partial-band 
interference,  multitone  jamming,  and  multiple-access  interfer¬ 
ence.  The  robustness  is  due  to  the  iterative  turbo  decoding  and 
demodulation,  the  channel  estimator  based  on  the  expectation- 
maximization  algorithm,  and  the  spectrally  compact  modula¬ 
tion. 
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